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Abstract. We study the prospects for detecting gamma-rays from decaying Dark 
Matter (DM), focusing in particular on gravitino DM in R-parity breaking vacua. 
Given the substantially different angular distribution of the predicted gamma-ray 
r*j \ signal with respect to the case of annihilating DM, and the relatively poor (of order 

Qh' 0.1°) angular resolution of gamma-ray detectors, the best strategy for detection is in 

q , this case to look for an exotic contribution to the gamma-ray flux at high galactic 

latitudes, where the decaying DM contribution would resemble an astrophysical extra- 
galactic component, similar to the one inferred by EGRET observations. Upcoming 
experiments such as GLAST and AMS-02 may identify this exotic contribution and 



discriminate it from astrophysical sources, or place significant constraints on the mass 
and lifetime of DM particles. 
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1. Introduction 
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A tremendous theoretical and experimental effort is in progress to clarify the nature of 
the elusive Dark Matter that appears to dominate the matter density of the Universe [TJ 
12]. The most studied DM candidates are Weakly Interacting Massive Particles, that 
achieve the appropriate relic density by freezing- out of thermal equilibrium when their 
self-annihilation rate becomes smaller than the expansion rate of the Universe. The 
characteristic mass of these particles is (9(100) GeV, and the most representative 
and commonly discussed candidates in this class of models are the supersymmetric 
neutralino, and the B*- 1 -* particle, first excitation of the hypercharge gauge boson, in 
theories with Universal Extra Dimensions. These particle are and will be searched for 
via collider, direct and indirect searches. In particular, the latter are based on the very 
same mechanism that controls the relic density of DM, i.e., self-annihilations. In fact, 
although the annihilation rate in the local Universe is on average severely suppressed, 
it can still be extremely high at the centre of dense DM halos, since it is proportional 
to the square of the DM particles number density. The prospects for indirect detection 
of annihilating DM have been extensively discussed (see [H 12] and references therein). 

However, self-annihilating relics are not the only DM candidates, and indirect DM 
searches are not only relevant for self-annihilating particles. Three of us have recently 
studied an excellent DM candidate, the gravitino in R-parity breaking vacua, that can 
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achieve the appropriate relic density through the thermal production in the early high- 
temperature phase of the Universe, and that naturally leads to a cosmological history 
consistent with thermal leptogenesis and primordial nucleosynthesis [3]. 

Since R-parity is broken, gravitinos can decay into a photon and a neutrino [4|, 
although with a lifetime that, being suppressed both by the Planck mass and by the small 
R-parity breaking parameters, is naturally much longer than the age of the Universe [5]. 
Similarly, for sufficiently small R-parity breaking, also neutralinos [SI U\ and axinos [8] 
are dark matter candidates which can decay into a photon and a neutrino. These 
scenarios thus predict a diffuse flux of photons and neutrinos that, by comparison with 
existing observational data, can be used to set constraints on the mass and lifetime of 
the decaying particles. Interestingly, an excess in the galactic component of the diffuse 
gamma-ray flux measured by EGRET has been claimed in Ref . [9] , at energies between 
1 and 10 GeV. A more careful analysis of the Galactic foreground has led Strong et 
al. to a new estimate of the extra-galactic component [10] with a significantly different 
spectrum with respect to the previous analysis. More recently, Stecker et al. [llj pointed 
out a possible error in the energy calibration of EGRET above 1 GeV, a circumstance 
that if confirmed would make any interpretation of EGRET data in terms of exotic 
components, such as DM annihilation or decay, unreliable, if not meaningless. 

In view of these and other systematic uncertainties [12], we will not try here to 
fit the EGRET data with the gamma-ray flux produced by decaying DM, although 
we regard this coincidence as interesting and deserving further attention. We perform 
instead a careful analysis of the signal that might be detected with the next generation of 
gamma-ray experiments. Similar analyses has previously been carried out for decaying 
DM candidates with masses in the keV range, such as a scalar modulus [13J or a sterile 
neutrino [TH [T5] . More recently also the case of small mass splittings and heavy DM 
decaying into MeV photons has been discussed in order to explain the COMPTEL excess 
in the photon flux [TSJ E] . 

Unlike the case of stable neutralinos and other WIMPs, the rate at which gravitinos 
produce photons is proportional to the density of DM particles, as appropriate for 
decaying DM particles, not to the square of the DM density. As a consequence, the 
strategies for indirect detection must keep into account the different angular distribution 
of the predicted signal, and the different ratio between galactic and extra-galactic 
contributions. Although the situation is very similar to the case of other decaying 
DM candidates, such as decaying sterile neutrinos, the angular resolution of experiments 
sensitive to photons from decaying gravitinos, typically above 5 GeV [3] , are much worse 
than X-ray telescopes, relevant for sterile neutrinos. Here we study the best strategies 
to detect an exotic component in the gamma-ray diffuse flux with future experiments 
such as the upcoming gamma-ray satellite GLAST, scheduled for launch in the next few 
months, and with AMS-02. Although in our analysis we adopt gravitinos as our fiducial 
DM candidates, our results can be applied to any decaying DM particle in a similar 
range of masses into monochromatic photons. 

The paper is organised as follows: in Sec. 2 we discuss our expectation for 
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gravitino lifetime and the decay channel, while in Sec. 3 the Galactic and extra-galactic 
contributions to the gamma-ray flux from gravitino decays. In Sec. 4 we consider 
alternative targets for indirect detection such as dwarf galaxies and galaxy clusters. In 
Sec. 5 we discuss our results and compare the indirect detection strategies of decaying 
and self- annihilating DM. Finally we give our conclusions in Sec. 6. 



2. Gravitino decay 

As is the case for the proton, we do not know if the DM particle is absolutely stable. In 
the case of supersymmetric candidates, usually R-parity is invoked to make the proton 
sufficiently long-lived and it automatically gives that the Lightest Supersymmetric 
Particle (LSP) is stable. On the other hand R-parity is not the only possibility for 
protecting the proton from rapid decay and a small amount of R-parity violation does 
not rule out the possibility of having supersymmetric DM if the LSP is very weakly 
interacting with the R-parity violating sector or the decay is highly suppressed by 
phase space. As an example, in [3], we have considered a supersymmetric extension 
of the standard model with small R-parity and lepton number violating couplings and 
a gravitino LSP. The model predicts a small photino-neutrino mixing \U^ V \ = O(10~ 8 ), 
which leads to the decay of the gravitino into photon and neutrino [5], 

Using Mp = 2.4 x 10 18 GeV, one obtains for the gravitino lifetime 

2-body o Q lfl27 (\U^\\ 2 ( ?"3/2 \~ 3 ( ] 

r 3/2 _ 3.8x10 { 1QGeV ) ■ ( 2 ) 

At tree level this decay channel can be suppressed if the sneutrino v.e.v. responsible 
for the photino-neutrino mixing is very small, but even if the mixing vanishes, the 
decay can take place via one loop diagrams. The loop induced decay has been recently 
computed in [18] , where it has been shown that also in this case the channel dominates 
over the 3-body decay into fermions [19] for small gravitino masses. For this reason 
we will concentrate on this particular channel and will assume in the following that our 
Dark Matter candidate decays into a photon and neutrino producing two monochromatic 
lines at energy equal to ttidm/^ with a lifetime of the order of 10 27 s or larger. On the 
other hand, if the gravitino is sufficiently heavy it could decay into W or Z bosons, 
producing through fragmentation a continuous spectrum of photons with a characteristic 
shape [20]. The neutrino flux in the few GeVs energy range is unfortunately overwhelmed 
by the atmospheric neutrino background and so its detection seems much more difficult 
than that of the gamma-ray flux. 

Note that the signal in gammas would be the same, only twice as strong, for the 
case of a scalar DM candidate decaying into two photons. 
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3. Gamma-Rays from DM decay 

If the DM particles decay all around us, we expect two sources for a diffuse background. 
We have the DM decaying in the Milky Way halo nearby and in addition those decaying 
at cosmological distances. 

Let us first consider the latter ones, which have been more intensively studied 
in the literature [5j [21] . The decay of DM into photon and neutrino at cosmological 
distances gives rise to a perfectly isotropic extragalactic diffuse gamma-ray flux with a 
characteristic energy spectrum, corresponding to a red-shifted monochromatic line. A 
photon with measured energy E = 1113/2/ (2(1 + z)) has been emitted at the comoving 
distance x( z )i with dx/dz = (1 + z)~ 3 / 2 /(a H JQM (1 + «(1 + z )~ 3 ))- Here a and H 



are the present scale factor and Hubble parameter, respectively, and k = Q\/Qm — 3, 
with + Qm — 1, assuming a flat universe. Then we obtain for the photon flux 

dh m DM V \m DM / I \ m DM ' V "\dm/ 



with 



a <> - - ^ m-T (sr • « 



4nT DM m DM H {l M 

here tdm is the DM particle lifetime and is given by Eq. (j2J) for the gravitino case. 
We have taken the particle density to be equal to the Cold Dark Matter density as 
^DMh 2 = 0.1, and the other constants as p c = 1.05 h 2 x 10~ 5 GeVcm -3 , total matter 
density Vt M = 0.25 and H = h 100 km s" 1 Mpc" 1 with h = 0.73 [22]. We are 
considering here tdm 3> Hq 1 so that we can neglect in the above formula the depletion 
of the number density due to the decay. 

In addition to the extragalactic signal there is an anisotropic sharp line from the 
halo of our galaxy with an intensity comparable to the extragalactic signal [13]. The 
flux from the decay of halo DM particles is given by the density profile, i.e. 
dJhaio . 2 / 2E \ 

—ttt = A halo 6 1 , (5) 

dE m DM \ m DM J 

where 

1 f -> 

Ahalo = ~. / Phalo(l)dl . (6) 

4:TTT DM m DM Jl.o.s. 

The ratio Ah a io/A eg only depends on cosmological parameters and the halo dark matter 
density integrated along the line of sight 0. Hence, the intensity and angular distribution 
of the halo signal is very sensitive to the distribution of the dark matter in the Milky 
Way. Surprisingly, for typical halo models, this ratio is of order unity [T3j . 
Consider a Navarro-Frenk- White profile for the DM matter of our galaxy, 

PNFw( r ) = —, — TT^T — / — ^ (7) 

\ The coefficients A^aio and A eg are related to the coefficients C 7 and D 7 in our previous paper [5] 
by a factor 2/7713/2- Let us also note that in that reference there is a typo in the definition of C 7 : it 
should say 10~ 7 instead of 10~ 6 . 
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Figure 1. Left: EGRET diffuse emission in the energy range E=[4,10] GeV. Right: 
Sum of the Galactic plus extra-galactic contributions to the gamma-ray ffux from 
gravitino decay for t dm = 4 x 10 27 s and thdm = 10 GeV and the NFW profile given 
in Eq. ©. 



with ph = 0.33 GeV/cm 3 = 0.6 x 10 5 p c giving the halo density normalisation and r c = 20 
kpc the critical radius where the profile slope changes. For any given point along the 
line of sight, the distance r from the centre of the galaxy can be expressed as function 
of the galactic coordinates, the longitude I and latitude b, and the distance from the 
Sun s in units of R & as 



r 



2f s,b,l) = R 2 Q (s - cos 6 cos/) 2 + (1 - cos 2 6 cos 2 /) . (8) 



The flux factor then reads 

Rq f°° 

A halo (bJ) = : / ds PNFw(r(s,b,l)) . (9) 

This expression can be used to give the flux dependence on the angle, as long as we 
are far away from the central cusp. In that region a more appropriate quantity is the 
average flux on the solid angle corresponding to the detector resolution around the 
direction (6, /), i.e. 

(A halo (b, = -L / alQ A halo (b, /) , (10) 
AS 2 J An 

where the infinitesimal solid angle is given by dVL = dlcos(b)db. We consider therefore 
the average flux for an angular resolution of 1° as measured by EGRET. The result of 
the numerical integration for the halo flux plus the isotropic extragalactic component 
is shown in the right panel of Fig. [Tj which illustrates the decrease of photon flux away 
from the galactic centre, both in longitude and latitude. Note that the dependence 
is not very strong and indeed Ahaio changes only by a factor 20 between the galactic 
centre and the anti-centre, and within a factor of 8 if one cuts out the region within 10 
degrees around the galactic plane. Averaging over all sky excluding the galactic plane, 
we obtain Ah a io/A eg ~ 0.76, so that the line is actually dominating the signal. In fact 



Gamma-Rays from Decaying Dark Matter 



6 



the total flux of diffuse gamma-rays for all sky directions is given by 

•«,,= /<mjf rf *(^ + ^(M)) (ID 

= 4nA eg — r— I ™s(b)db dlA halo {b,l) (12) 

Oy K J— 7r/2 J— 7r 

.4 



4nA eg fo.5 + (13) 



leg 



= 1.5 x 10- (cm' .)"' , (14) 

where we have used k = 3 and Eqs. (JH) and (jBJ). 

The expected signal from Dark Matter decay in the halo can be compared with 
the diffuse gamma-ray flux observed by EGRET. Contour lines of constant flux, with 
photon energies between 4 GeV and 10 GeV are shown in the right panel of Fig. [TJ As 
expected, the smallest flux is observed in the directions of the north and south poles, 
i.e. orthogonal to the galactic disk. According to Fig. [TJ the signal in these directions 
is larger than in the direction opposite to the galactic centre. We therefore conclude 
that the signal from the MW halo could be most effectively observed looking away from 
the galactic disk, which generates most of the background, in direction of the poles, in 
contrast with the strategy usually adopted for the detection of the self-annihilating DM 
signal. 

The photon spectrum is dominated by the sharp line coming from our local halo, 
while the red-shifted extragalactic signal is appreciably lower. The position of the 
line allows a direct measurement of the DM particle mass and the height is inversely 
proportional to the lifetime. In Fig. [2j we show the expected signal for decaying particle 
mass of 10 GeV and for a lifetime of 10 27 s in comparison to the extragalactic EGRET 
data [9J [10]. We mimic the finite energy resolution of the detector by convolving the 
signal with a Gaussian distribution and average the halo signal over a cone of 80° 
around the poles. The height and width of the line depend as usual on the energy 
resolution of the detector; here we have taken 15% as energy resolution, as quoted 
by EGRET in this energy range. Note also that the DM signal peaked at 5 GeV 
corresponds to the expectation of the model of gravitino LSP with R-parity and B-L 
breaking discussed in [3]. A word of caution though is in order in the comparison 
between data and signal: the EGRET extragalactic background displayed here has been 
extracted assuming isotropy, while our halo emission is mildly anisotropic away from 
the galactic plane. GLAST is expected to provide much better data at these energies, 
allowing a much more detailed analysis of the angular and spectral properties of the 
diffuse gamma-ray flux. 



4. Alternative targets 



The Milky Way has dwarf galaxies as satellites, which have a large mass to light ratio, 
like Draco and Ursa Minor. One may therefore hope that the flux from decaying DM 
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Figure 2. Energy spectrum of extragalactic and halo signal compared to the EGRET 
data. The data points are the EGRET extragalactic background as extracted by Strong 
et al. in [10] , while the short-dashed (blue) line shows the powerlaw fit from Eq. (|27|) 
obtained previously by [9]. The extragalactic and halo signals for tdm — 10 27 s and 
rriDM = 10 GeV are respectively the long-dashed (green) and dotted (magenta) lines, 
while the solid (thick red) line shows the sum of these contributions with a powerlaw 
background (thin red line), which has been obtained fitting the low energy EGRET 
points. 



is significantly enhanced in these directions in the sky like in the case of annihilating 
DM |23j. In the following we shall study the dependence of the enhancement on the 
angular resolution of the detector and the mass of the dark matter constituents. For 
simplicity, we use the isothermal profile, 

Phalo(r) = 1+ P r ° 2/r 2 , (15) 

for which one can easily derive simple analytic expressions for the photon flux. 

Integrating along the direction of sight and using Eqs. ([5]), §6§ and ([8]), one finds 
for the photon flux from decaying dark matter in the Milky Way halo 

1 1 p r 2 c 



Jhalo(b, I) 



47T T DM m DM Rq h _ cog 2 I CQS 2 I + r 2J R 2^ 








, | , cos b cos I , . 

x [ - + tan" 1 ( ; ) ) . (16) 

- cos 2 b cos 2 1 + r 2 / Rq t 

With Rq ~ 8.5 kpc, r c ~ 3.5 kpc and po — 1-37 GeV cm -3 , this yields in the direction 
of Draco (b D = 34°, l D = 86°; cos/ D cos6 D ~ 0.06): 

■WWiO - 0.8 x 10-' (em's *p)- (^) _1 ( I ^ F )"' • (17) 
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The photon flux observed from a distant dwarf galaxy crucially depends on the 
angular resolution of the detector. Averaging over a cone with small opening angle S, 
directed toward the centre of the dwarf galaxy, one has 

1 1 f S rr m ax(a) I 

2T DM m DM AiUv Jr min (a) r z + a z r z 

here AQ = n5 2 is the infinitesimal solid angle, r(a) = JJd — r) 2 + a 2 r 2 , d is the distance 
to the dwarf galaxy, p^g is the isothermal profile in Eq. (1151) ; we have also taken a finite 
tidal radius r m into account, which leads to the finite integration domain given by 
t{ci) 2 < r^. Performing the integrations, one obtains in the relevant case where 5< 1 
and Sd, r c <C r m , 

Jd ( 5 ) ~ Porc ( 2rc 2rc | ) , (19) 

9 ^T DM m DM \r c + Jr 2 + 5 2 d 2 nr m 



where terms 0(5 2 ) and 0(8dr c /r^ n ) have been neglected. Up to corrections 0(r c /d), 
the numerator of the prefactor is precisely the line of sight integral of the dark matter 
profile for r m — > oo, 



TWc = / drp dg (r) , (20) 

J l.o.s 

and for 5, l/r m — > 0, the bracket becomes one. 

Taking as an example Draco, typical parameters are r c = 0.1 kpc, p$ r = 
28.4 GeV cm" 3 , r® r = 1.7 kpc, d Dr = 80 kpc [15J. The correction due to the finite tidal 
radius in Eq. ( |T9i) is then negligible, and in the line of sight approximation (5 = 0) one 
obtains the flux 

M0) <>> 2 x 10- 7 (cm 2 s str)" 1 ( -^V' ( -^Q" 1 , (21) 



V10 27 s/ VlOGeV, 

which is about three times larger than the flux from the halo in the same direction, as 
pointed out in [15]. The correction factor for finite opening angle in Eq. (|T9|) can only 
be neglected for 5 < r c /d = 1.3 x 10 -3 , which corresponds to the angular resolution 
of GLAST of about 0.1°. For 5 ~ 0.004 the signals from Draco and the halo have 
equal strength, but the corresponding field of view is AQ ~ 5 x 10 -5 , yielding the 
flux $ D (0.004) ~ 3.6 x lO-^cn^s)- 1 , too small to be observed even by GLAST [Mj. 
We conclude that for dark matter particles with masses in the GeV range the flux 
enhancement in the direction of dwarf galaxies is currently not of interest. This is 
different for masses in the keV range (15] where the flux is six orders of magnitude 
larger for the same DM density and in which case the angular resolutions of the X-rays 
detector is much better than for gamma rays. 

Another potentially important source of gamma rays is the Andromeda Galaxy, due 
to its proximity, large apparent size and privileged position in the sky away from the 
Milky Way centre. To estimate the photon flux received from the Andromeda Galaxy, 
we approximate the dark matter distribution by the isothermal profile, Eq. ( jT5l) . with 
r^ 31 = 1.5 kpc, p™ 31 ~ 15.7 GeV cm" 3 , < 31 = 117 kpc [25]. The total flux received 
from a cone directed toward the centre of the Andromeda Galaxy with a small opening 
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angle 5 is given by Eq. (1191) . On the other hand, for the angular resolution of GLAST, 
the gamma ray flux from the Milky Way halo in the direction to the Andromeda Galaxy 
(&M3i = — 22°,/m3i = 121°) is approximately 

Jwo(&M3i,/M3i) * 0.5 x 10- 7 (cm 2 s str)" 1 (^)~ ' (t^^) 1 ■ ( 22 ) 

The opening angle at which the Milky Way signal equals the Andromeda Galaxy signal 
is around 5°, which corresponds to a flux $M3i(5°) — 1.5 x 10~ 9 (cm 2 s) _1 . This is above 
the sensitivity of detection of GLAST (although below the EGRET sensitivity), hence 
we conclude that GLAST might be able to see the Andromeda Galaxy as a gamma 
ray source over the background from the decaying dark matter in the Milky Way halo. 
However, one should also note that the extraction of an extragalactic gamma ray flux at 
such low galactic latitudes is intricate, and the Andromeda signal from dark matter decay 
could be easily masked by photons from standard astrophysical processes occurring in 
the Milky Way disk. 

Finally, let us discuss other potentially interesting sources of gamma rays, namely 
nearby galaxy clusters located at high galactic latitudes, like Coma or Virgo. The 
photon flux can be computed along the same lines as for dwarf galaxies, where the 
density profile Eq. ( fl5i) has to be replaced by the profile for the isothermal /3-model 
[261 [271 

d _ d 3 + r 2 /r 2 

P ~ Po (1 + r2/r 2 )2 , (2d) 

that approximately describes the distribution of dark matter in a galaxy cluster. The 
result for 5 <C 1 and 8d, r c <C r m can be approximated by: 

U*)* ^ ( i ^ - — + •••! • (24) 
2r DM m DM yJrl + Pd 2 7rr m J 

Nearby galaxy clusters have a large angular size and could not appear as point 
sources. To be precise, the core of the Coma Cluster has a size of 0.3 Mpc compared to 
a distance of 98 Mpc, corresponding to 0.17° of angular size, which is larger than the 
angular resolution of GLAST. On the other hand, the Virgo Cluster has a core radius 
of 10 kpc and lies at 18 Mpc, which translate into an angular size of 0.03°, slightly 
smaller than the angular resolution of GLAST. In consequence, for these objects the 
line of sight approximation 5 = is not a good approximation and in order to obtain a 
reliable estimate for the photon flux the complete expression Eq. f[2"4"l) has to be used. 
Taking p$ = 10~ 2 GeV cm -3 and p^ = 2.3 GeV cm -3 for the Coma and Virgo clusters, 
respectively, we obtain for the angular resolution of GLAST 

J Coma (0.1°) ~ 4 x 10- 7 (cm 2 s str)- 1 < T ° M ^ < ™ DM ' " 



10 27 sj V 10 GeV 
i / ™ \ -l 



W0.1-) = 10-(cmW)->(^)-(^)- . (25) 
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which are around one order of magnitude larger than the photon flux from the Milky 
Way halo in the direction of the galactic pole, 

J halo (b = tt/2) - 0.7 x 10- (cm 2 s str)" 1 (^f-)" ■ (26) 

Nevertheless, the total flux received from these objects is so small that in order to 
distinguish a gamma ray signal from the Coma and Virgo clusters, it would be necessary 
a sensitivity of 7 x 10~ n photons cm~ 2 s _1 , which is more than one order of magnitude 
lower than the GLAST sensitivity [24] . 

5. Discussion 

Should the gamma-ray line from DM decay be discovered at high galactic latitude, the 
problem will arise of how to discriminate it from scenarios with self-annihilating DM. 
The presence of a spectral line in the diffuse signal at high galactic latitude has in fact 
been proposed, several years ago, as a signature of annihilating DM particles [28] (see 
also Refs. [291 ED]). We note however that there are at least four important differences 
between the two scenarios: 

• the angular profile of the predicted gamma-ray signal; 

• the comparison between galactic and extra- galactic components; 

• the ratio between the line and the gamma-ray continuum; 

• the angular power spectrum of the gamma-ray background. 

We stress in fact that in the case of self-annihilating DM, the gamma-ray flux from 
a given direction in the sky is a steeply falling function of the angle with respect to the 
Galactic centre, and it is in particular much steeper than the flux from decaying DM, as 
can be seen in Fig. El That is precisely why we have focused on high galactic latitudes 
in Sec. 2. Furthermore, the extra-galactic component of the gamma-ray flux in the 
case of annihilating DM is unlikely to be detected in absence of a strong signal from 
the Galactic centre [31] . In other words, in order for an extra-galactic component to be 
detected at high latitudes, the gamma-ray signal from the galactic centre should be easily 
detectable, despite the much stronger astrophysical background toward the innermost 
regions of the Galaxy. Thus, in case a line is observed at high galactic latitude, the 
absence of a similar line in the gamma-ray spectrum of the Galactic centre would favour 
an interpretation in terms of decaying DM. 

It is still possible that astrophysical processes, such as the formation of spikes [321 
[331 [34] and mini-spikes [351 ESI EZ], i.e. large DM overdensities around Supermassive 
and Intermediate Mass Black Holes respectively, modify this picture, by boosting 
the gamma-ray signal from cosmological structures much more than the flux from 
the Galactic centre [381 139] - Even in this case, however, it should be possible to 
discriminate between the two scenarios, by studying the spectral features of the extra- 
galactic background. In fact, for most self- annihilating DM candidates, such as the 
super symmetric neutralino and the particle in theories with Universal Extra 
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Figure 3. Angular profile of the gamma-ray signal as function of the angle 9 to the 
centre of the galaxy for a NFW halo distribution for decaying DM, solid (red) line, 
compared to the case of self-annihilating DM, dashed (blue) line. Both signals have 
been normalised to their values at the galactic poles, 9 — ±90°. The central cusp 
is regularised by assuming in both cases the GLAST angular resolution of 0.1° and 
integrating on the solid angle as in Eq. (flT)]) . 

Dimensions, direct annihilation to photons is severely suppressed with respect to other 
channels such as annihilation to quarks or gauge bosons. It follows that the annihilation 
spectra are characterised by a continuum emission that is inevitably associated with 
the line signal, as shown in Fig. HI Furthermore, masses below 50 GeV are usually 
considered unlikely for annihilating candidates. For instance, the current constraint in 
the neutralino mass is 40 GeV (assuming unification of gaugino mass parameters at 
the GUT scale) [22] while electroweak precision data exclude B^ 1 ) masses below 300 
GeV [10]. 

Although the relative importance of the line can be in some cases particularly 
high, as e.g. in the case of Inert Higgs DM |41j . the peculiar shape of the angular 
power spectrum of the gamma-ray background |42j can be used as a diagnostic tool 
to discriminate between annihilating and decaying DM. As for all the other strategies, 
this method heavily relies on the assumption that a sufficient number of photons are 
collected above the astrophysical background to allow a statistically meaningful analysis. 
In the worst case scenario, null GLAST searches can be used to exclude regions of the 
parameter space that lead to observable fluxes. 

In Fig. [5] we show an exclusion plot in the (m DM ,T DM ) plane, where we show the 
regions of the parameter space that already are ruled out by a comparison with EGRET 
data. In order to investigate the discovery potential of GLAST, we refer to the official 
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Figure 4. Spectrum of decaying DM extra-galactic component, solid (red) lines, 
compared with the spectrum of annihilating DM (continuum as long-dashed (green) 
plus line in short-dashed (blue)). For the case of decaying DM, E max = mjjjif/2, 
whereas for the case of annihilating DM, E max — ttidm- 



publications of the collaboration [24J where the sensitivity to astrophysical lines as a 
function of energy is studied for two different sky regions, namely an 'annulus' region 
where the signal to background ratio is maximised for annihilating DM candidates, and 
a 'high latitude' region, defined as the region of the sky with galactic latitude higher 
than 20 degrees and excluding a 35 degrees circle around the Galactic centre. In our 
case, there is a considerable reduction of the background, and a not-so-large decrease 
of the signal when going at high galactic latitude, for the reasons discussed above. So 
it is not a surprise that the 'high latitude' region is more favourable for constraining 
decaying DM. In Fig. Owe show the reach of GLAST for the two regions along with the 
EGRET-excluded region. For the present bound we take conservatively the requirement 
that the peak in the energy spectrum with 15% energy resolution remains below the 2a 
band of the EGRET spectrum obtained by Sreekumar et al. in [9], 

j j / TP \ -2.1±0.03 

dJ egret = (7.32^0,34) xl0- 6 (cm 2 sstr GeV)- 1 ( =— -) .(27) 

dE v ; v ; V0.451 GeV/ K ' 

The different mass dependence of the EGRET bound is due to the fact that it is a 

constraint on 4g instead than on the integrated flux and on the energy dependence of 

the GLAST sensitivity. The limit disappears at 100 — 120 GeV, which is the maximal 

energy plotted in the EGRET fit, but note that data above 10 GeV have large errors 

and were not used to obtain Eq. (!2Tj) . 
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m DM [GeV] 

Figure 5. Exclusion plot in the (thdm ,t~dm) plane based on EGRET data. All the 
region below the solid (red) line is excluded by the requirement that the halo line is 
below the measured flux as explained in the text. We also show the region of the 
parameter space where GLAST could discover the annihilation line: the short-dashed 
(blue) line corresponds to the GLAST sensitivity in the 'annulus' region, while the 
long-dashed (green) line to the 'high latitude' region. 

6. Conclusions 

We have discussed the prospects for indirect detection of gravitino DM in R-parity 
breaking vacua with upcoming experiments such as GLAST and AMS-02. The search 
strategy is in this case significantly different with respect to annihilating DM particles, 
due to different angular profile of the predicted gamma-ray signal and the different 
relative importance of the extra-galactic background. We found that the predicted 
signal from DM decays in the galactic halo would resemble an extra-galactic component 
with a shape and normalisation similar to the one inferred by EGRET data. 

Without trying to fit the controversial EGRET data with our model, we have 
determined the regions of the (mM^M) plane where GLAST could discover the 
gamma-ray line from gravitino decay, and we discussed how to discriminate this signal 
from astrophysical sources and from a signal originating from annihilating DM. This 
discrimination is based on the fact that: the line may appear at energies below 50 
GeV, which is somewhat challenging for popular annihilating DM candidates; it would 
exhibit no continuum flux (unless the mass of the gravitino is above the W mass); there 
would be a weak line signal from the galactic centre, possibly hidden in the galactic 
background; the angular power spectrum of the signal would be different from the case 
of annihilating DM, and from extra-galactic astrophysical sources. 
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In the worst case scenario, i.e. in case of null GLAST searches, the data can be 
used to set constraints in the (m DM ,r DM ) plane and improve the bound on the DM 
lifetime into photons by more than one order of magnitude, actually even more than 
two for masses below 1 GeV. 
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